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Abstract: This study presents a tomographic survey of a subset of the outer halo (10-40 kpc) drawn 
from the Sloan Digital Sky Survey Data Release 6. Halo substructure on spatial scales of > 3 degrees 
is revealed as an excess in the local density of sub-giant stars. With an appropriate assumption of a 
model stellar isochrone it is possible for us to then derive distances to the sub-giant population. We 
describe three new candidate halo substructures; the 160- and 180-degree over-densities (at distances 
of 17 and 19 kpc respectively and radii of 1.3 and 1.5 kpc respectively) and an extended feature at 28 
kpc that covers at least 162 square degrees, the Virgo Equatorial Stream. In addition, we recover the 
Sagittarius dwarf galaxy (Sgr) leading arm material and the Virgo Over-density. 

The derived distances, together with the number of sub-giant stars associated with each substructure, 
enables us to derive the integrated luminosity for the features. The tenuous, low surface brightness 
of the features strongly suggests an origin from the tidal disruption of an accreted galaxy or galaxies. 
Given the dominance of the tidal debris of Sgr in this region of the sky we investigate if our observations 
can be accommodated by tidal disruption models for Sgr. The clear discordance between observations 
and model predictions for known Sgr features means it is difficult to tell unambiguously if the new 
substructures are related to Sgr or not. Radial velocities in the stellar over-densities will be critical in 
establishing their origins. 
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1 Introduction 

How does a galaxy such as our Milky Way form? This 
remains a pressing question for moder n astrophysics. 
The seminal work of Eggen et al. ( 1962 ) proposed that 



galaxy formation occurred in the wake of a monolithic 
collapse of an isolated proto-galactic cloud. Searle <£] 
Zinn| ( 1978 ) challenged this picture of isolated forma- 



tion by proposing instead that the halo was assembled 
by the accretion of numerous small entities over an 
extended period. 

Recent observational evidence suggests both for- 
mation scenarios have a role to play. | Car olio et al.| 



(2007) presents evidence that the inner halo (that of 



galactocentric radius < 15kpc) is dominated by highly 
eccentric, prograde orbits and a metallicity of [Fe/H] 
^ —1.6 due to in-situ formation, whereas the outer 
halo shows a more uniform distribution of eccentrici- 
ties, includes highly retrograde orbits, and lower mean 
metallicity ([Fe/H] ^ —2.2) that, they argu e, has an 
accretion origin. The [Kinman et al. ( 2007 ) study of 
the local halo, similarly identifies a population that 
possesses retrograde rotation and streaming motions 
(i.e. low velocity dispersion) and another population 
that exhibi ts negligible rotation and less indication of 
streaming. Kinman et al. (2007| show that the hor- 



izontal branch morphologies of these two halo com- 
ponents are consistent with those seen in the young 
and old globular cluster systems respectively. This re- 
inforces the conclusion of |Zinn| (1985 ) that the inner 
halo is consistent with formation from rapid collapse 
and that the outer halo was formed from ongoing ac- 



cretion of dwarf galaxies. 

Simulations of galaxy formation in the context of 
the prevailing ACDM cosmology predict that the ma- 
jority of stars in the outer halo formed in progenitor 
dwarf gal axies that were subsequently accreted (s ee for 
example, [Freeman Bland-Hawthorn (2002)). |Bul- 



lock &: Johnson| ( |2005| ) show that spatially coherent 
substructures formed of tidal debris should survive in 
the outer halo for several Gyrs and that their kine- 
matic signatures (in velocity space) should remain de- 
tectab le well after spatial substructure is no longer ev- 
ident (Helmi et al. [2003 ). On the other hand, a halo 
that formed from in-situ star formation would possess, 
at the current epoch many dynam i cal tim es later, neg- 
ligible substructure. Bell et aT] (pOOSl compare the 
level of substructure seen in the SDSS data with cos- 
mological simulations (over galactocentric distances 1- 
40 kpc) and find that the Milky Way's halo is consis- 
tent with assembly from dwarf galaxies. 

The most prominent contribution to outer halo 
substructure results from the on-going disrupt ion of 
the Sagittarius Dwarf Spheroidal galaxy (Sgr, |Ibata 



et al.|p!994 ). The leading and trailing arms of Sgr are 



2003 New- 



seen to encircle the sky ( Majewski et al 

berg et al.l[2002l [2007) [Belokurov et al.||2006| " |keller 
et al.||20'Q8 ). The extensive debris of Sgr offer an im- 



port ant prob e of the mass and shape of the Galaxy 
( Ibata et al.p OOl; Hel mi Wh ite 2001; Helmi 2 004| 
^ohnston et al.n2005, ,Law et al.|2005, ,Fellhauer et al.| 



2006) 



Another outer halo substructure, the Virgo Over- 
density (VOD) was found as an over-density in RR 
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Lyrae variables (RRLs, [Vivas et al.||2001| [2004| jvT] 
IKeller et al. 2008) and in SDSS as 



a diffuse structure of main-sequenc e stars ( .Juric et al.^ 



2008| [Newberg et~aL]|2002| |2007|) and blue horizon- 
tal branch stars (BHBs, Ivezic et al. |2005|). Radial 



velocity measurements by|Duffau et al.| ( |2006|) 



New- 



berg et al.| ( |2007| ), [Vivas et al.| ( |2008| ), and jPrior et al 



(|2009[) reveal that a subset of associated RRLs share a 
common spatial velocity. This moving group is termed 
the Virgo Stellar Stream (VSS) to distinguish it from 
the general stellar over-density and potentially other 
velocity groups that may spatially coincide. 

The study of|Duffa u et"aL ]([2006| explored the spa- 
tial extent of the VOD by considering the excess of the 
observed luminosity function over an off-source control 
field. This study discerned a size for the VOD of at 
least 106 s quare degrees. A s imilar technique is imple- 
mented by [Prior et al. (2009) to derive a spatial extent 
for the VQD of ^ 760 square degrees. The study by 
Juric et al. (2008| utilized the photometric parallaxes 



of F-type main-sequence stars to describe a large dif- 
fuse over-density covering ^ 1000 square degrees and 
spann i ng a range of distance from 6-20 kpc. [Vivas [ 



et al. (2008| find a small concentration of stars with 



similar radial velocities to the VSS reaching to 12 kpc 
from a previous detection at 19 kpc. 

Martmez-Delgado et alT] ( [2004[ [2007[ ) suggest that 



the VOD is the result of the confluence of the lead- 
ing an d trailing arms of Sg r. However, as pointed 
out by [Newberg et al.[ ( [2007t and [Vivas et al.[ ( [2008[ ), 
models of the Sgr debris stream predict highly nega- 
tive radial velocities contrary to the +100-130 kms~^ 
seen in the VSS. Extant models also fail to match the 
observed de nsity enhancement of t he VOD. N ewberg] 



et al. (2007) and Vivas et al. (2008) conclude that the 



VOD is most likely tidal debris from a separate and dis- 
tinct accretion event. Recently, [Casetti-Dines cu et al.[ 
(2009) have presented a preliminary orbit for the VSS 



based on one RRL that places the system near peri- 
centre on an eccentric and highly destructive orbit. 

In light of the considerable uncertainty regarding 
the apparent size, distance and origin of the VOD we 
have presented a series of new obs ervational studies 
targeting this region. K eller et al.[ (|2008 ) has exam- 
ined the spatial distribution of RRL candidates in the 
VOD region and finds in addition to the over-density 
seen by the QUEST team at 17 kpc, a second region 
of over-density 15 degrees to th e south-east and at a 
distance of 20 kpc. The study of [Prior et a"L| ( [2009[ ) ex- 
amines the kinematics of a sample of RRL candidates 
from Keller et al. (2008). In contrast to previous stud- 
ies, this study finds RRLs that possess radial velocities 
consistent with membership of Sgr trailing arm debris. 
They conclude that a significant contribution from Sgr 
can not be ruled out at the present time. [Keller et al.[ 
([2009 ) utilizes the luminosity function exc ess over that 



expected from the Besancon galaxy model { Rob in et al. 
[2003 J to map the extent of the VOD. This study finds 
that the region consists of three spatially distinct over- 
densities that extend from declinations +2° to —15°. 

In the present study, we utilize the sub-giant popu- 
lation (those stars between the main-sequence turn-off 
and the base of the red giant branch) to perform to- 



mography of the outer halo in order to describe the 
VOD in more detail and to look for additional sub- 
structure. In Section [2] we present our method for 
determining the significance of, and distance to, halo 
substructure. Section [3] presents our results on the ex- 
tent and luminosity of a series of halo sub-structures 
with comparison, where possible to existing literature. 
Finally, in Section [5] we discuss our map of substruc- 
ture in relation to the tidal debris of Sgr. 



2 Sub-giant tomography 

Our technique is to take an identifiable stellar popu- 
lation and use this population as a tracer of spatial 
density by translating apparent magnitude to distance 
for the population. As discussed previously, this tech- 
nique has utililsed RRLs, BHBs, M giants and main- 
sequence (MS) turn-off stars. Here we utilise the sub- 
giant population. 

There are a number of advantages to the use of 
sub-giants. Firstly, sub-giants are a numerous popu- 
lation with which to trace substructure. Commonly 
used distance indicators like RRLs, BHBs and M gi- 
ants while presenting a population free from high con- 
tamination, are too sparse (typical values are 0.1 per 
square degree for M-giants to 10 per square degree in 
the case of BHBs) to present substantial contrast over 
their surrounds to probe spatial structure on small an- 
gular scales. 

The MS population presents a much more n umer 
ous populat ion (^ 7500 per square degree; see Juric 



et al. 



(2008)). Plentiful as they are, deriving distances 



to MS stars is fraught wit h considerable difficulty. As 
discussed in Juric et al.[ ([2008 ) conversion from the 
apparent magnitude to the distance of a MS star re- 
quires a well described photometric parallax relation. 
The technique requires highly accurate colours in order 
to precisely derive an absolute magnitude. As [Juric[ 
et al. (2008 their figure 7) shows, distances to stars 
at the MS turn-off possess a fractional distance uncer- 
tainty of between ±15— ^ 65%. This is a result of the 
extremely steep absolute magnitude - colour relation 
inherent in the photometric parallax method. In ad- 
dition, considerable systematics enter with the choice 
of a photometric parallax relation (up to ±70% in dis- 
tance for the bluest halo stars). 

Sub-giants, on the other hand, are less numerous 
but are significantly cooler stars than those at the MS 
turn-off. Furthermore, these objects result from a evo- 
lutionary phase of limited duration, as opposed to the 
MS. We can therefore, more precisely derive their dis- 
tance by relating their apparent brightness to that pre- 
dicted by stellar evolution isochrones. 

Sub-giants do, however, suffer from an ambiguity 
since they are selected from a window of the colour- 
magnitude diagram (CMD) that is inhabited by core 
helium burning stars on the red clump and horizon- 
tal branch, and lower MS stars. This has the poten- 
tial to lead to three peaks in measured stellar den- 
sity as a function of magnitude from a population of 
set distance as each of three populations are encoun- 
tered at intervals of approximately 3 magnitudes. If 
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populations of similar colours to the sub-giant branch 
are present they will lead to erroneous conclusions as 
to the distance of such over-densities. A population 
who's horizontal branch and/or red clump appears in 
the sub-giant window will be assigned a distance mod- 
ulus 3 magnitudes, or 4 times too close. We bear this 
caveat in mind in the discussion to follow. 

The region we have chosen to investigate in de- 



16 



tail is the "F ield of Streams" described by |Belokurov| 



et al. (2006). This region extends from declinations 
+34" to"^ and covers Right Ascension 120° to 220°. 
The area encloses the northern leading arm of the Sgr 
Stream and the VOD. We examine the significance of 
sub-giant excess over the nominal field population in 
square ^pixels' drawn from the photo metry of SDSS 
DR6 dAdelman-McCarthy et al.|2008| ). 

2.1 The determination of the signif- 
icance of sub-giant excess 

Our procedure of determining the significance of sub- 
giant excess at a position on the sky is as follows. We 
determine the number of objects within a window on 
the colour-magnitude diagram (CMD) that follows the 
locus of the globular cluster M3 sub-giant population 
between 0.35 < (g — r)o < 0.60 and is 0.2 magnitudes 
in height as shown in Figure^ The box is then placed 
at go magnitudes from 18.0 to 21.7 in 0.1 magnitude 
steps. Fainter than ^o=21.7 photometric uncertainties 
scatter halo turn-off stars into our sub-giant selection 
window. This reduces the contrast of possible super- 
imposed substructures. The bright limit of our magni- 
tude range is set to avoid contribution for the thin and 
thick disk main-sequence stars to our sub-giant sam- 
ple. This defines the effective distance range for the 
present st udy o f 10-40 kpc. As seen from simulations 
of Rob in et al. (|2003) (and furthermore, by observa- 
tions by [Conn et al.| (2008)), halo stars are dominant 
for go > 18.0. We then subtract a low-order polyno- 
mial fit to the number of sub-giant stars as a function 
of go magnitude. This leaves us with departures from 
the general background as shown in Figure [2] 

The spatial pixel size is 3° x 3°. Such pixels con- 
tain typically 200 stars within the colour selection win- 
dow at each magnitude step (hence a Poisson uncer- 
tainty of < 7% per magnitude step). Making the pixel 
size smaller introduces excessive shot noise and making 
the pixel larger removes spatial information. The rel- 
atively coarse grid size leaves us sensitive to extended 
stream structures but unable to resolve "fine-grain" 
(i3° wide) structur es such as the Or phan Stream (Be-^ 
lokurov et al.|2007t or the streams of |Grillmair| ( |2009| ) . 

We then consider the region 180° <RA< 190° and 
24° <Dec< 34° that, as evidenced by the CMD of the 
region, does not contain significant over-density. As 
can be seen from Belokurov et al. (2006) the northern 
branch of the bifurcated Sgr leading arm does pass 
through this area. However, it is sufficiently distant 
that it does not significantly enter the distance range 
probed by our sample. The quoted distance to this 
Sgr debris is, we contend, underestimated due to an 
assumed high metallicity as discussed below. From this 
large area we draw multiple samples that match the 
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Figure 1: An illustration of our selection of sub- 
giant stars from the field containing the globular 
cluster M3. The selection box (0.35< {g — r)o < 
0.60 and 0.2 mag. in height) shown is placed be- 
tween 18. 0< go <21.7 mag. at 0.1 mag. steps and 
the number of stars so selected recorded as a func- 
tion of magnitude. 



sample size typical of our 3° x 3° fields. This enables us 
to quantify the stochastic variation inherent in Figure 
[2] We can then transform the measured excess of sub- 
giants over the smooth field to a significance of sub- 
giant excess as a function of magnitude. 

By recourse to theoretical isochrones, we can relate 
the go magnitude at which the excess is observed to a 
distance. We have assumed an isochrone of 12 Gyrs 
and [Fe/H]=-1.7 (Girardi et al . 2004) that is typical 
of the population ex pected from a dwarf spheroidal 
galaxy ( [Mateo 1998) and is motivated by our previ- 
ous metallicity determina tions for RR Lyr aes in the 
VOD (average of-1.7dex, [Prior et al. 1200 91). Further- 
more, this is in line with the recent determination of 
the metallicity in the vicinity of the VOD by [An et al.[ 
([2009) who find [Fe/H]=-2.0=b0.5dex. 

As shown in Figure [2] for M3, we derive a distance 
modulus of 15.10=b0.02. This is in good agreement 
with the distance modulus of 15.09 from Harris ( 1996|. 
The distance modulus derived is only slightly depen- 
dent on the choice of age (e.g. a choice of 10 Gyrs 
results in a distance modulus of 15.12) but moderately 
dependent on the choice of metallicity, for example, 
assuming [Fe/H]=-1.3 results in a distance modulus 
of 14.94; [Fe/H]=-0.7 results in a distance modulus of 
14.52. 
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of 20 kpc and extends to our survey limit of 40 kpc 
at RA^ 195°. We can compare our findings with the 
locii o f the leading arm presente d in Belokurov et al.| 



Figure 2: A plot of the number of excess sub-giant 
stars as a function of magnitude resulting from the 
selection process described in Figure[T] The promi- 
nent peak seen corresponds to selection of the M3 
sub-giant branch. Using Monte-Carlo simulations 
we are able to quantify the excess of sub-giant stars 
in terms of statistical significance (right-hand axis, 
see text for details). 



3 Sub-structure revealed 

The significance of sub-giant excess as a function of 
distance was computed in each 3° x 3° x 0.6kpc 'voxel' 
in our survey area. The result is a data cube of over- 
density significance. For the purposes of illustration 
we have projected the data cube along the shortest 
axis (Declination). To form this projection we take 
the sum of the significance along the 12 pixels of the 
Declination axis. 

A number of sources were seen that correspond to 
single spatial (i.e. RA,Dec) pixels and confined to a 
narrow range in distance. These are due to known 
globular clusters in the field. In Figure [3] we show the 
the CMDs for four fields containing prominent globular 
clusters. The distances to three of the globular clus- 
ters are in good agreement with those in the literature 
(10.4, 14.4, 16.5 kpc compared with 10.4, 15.9, IT.lkpc 
from Harris|199 6, for M3, NGC5466 and NGC5024+5053 
respectively) . The CMD for Pal3 appears at a distance 
of 24 kpc, a factor of fo ur closer than its literature value 
(92.7kpc, Harris|p"996 ) due to the inclusion of its red 
clump in the sub-giant selection window. 

In the following analysis of the data-cube collapsed 
along the Dec axis we have excluded the RA, Dec pix- 
els that contain galactic globular clusters (this excises 
1.7% of the survey volume and hence is unlikely to 
modify our conclusions). Figure [4] is dominated by the 
Sgr leading arm. We utilise this feature as a further 
useful check of our methodology. The leading arm en- 
ters our survey volume at RA=122° and at a distance 



(2006) and Newberg et al. (2007) 



The distances in the study of Belokurov et al. ( 2006 ) 



are calculated from the relative offset between the colour- 
magnitude diagram (CMD) of the debris and that of 
the main body of Sgr. Since the most prominent fea- 
t ure of the CMD is the s ubgiant branch, the technique 
of Belokurov et al. ( 2006j) is analogus to that used here. 
However, the main body of Sgr has a mean metallic- 
ity of [Fe/H]=-0.4 ([Smecker-Hane McWiniam"[2Q02 



[Monaco enal|2005| |Chou et al.||2007[ with a con 
siderable low metallicity tail to Fe/H ~ —1.6) and 
the debris has mean [Fe/H] =-1.1 (but with reduced 
spread in Fe/ H relative to the core , |Chou et al.|2007t . 
Isochrones of Girardi et al. ( |2004 ) show that this dif- 
ference in mean metallicity leads to potentially a 0.4 
mag. under-estimation in the derived distance modu- 
lus of the Sgr leading arm material. If we apply this 
shift to the BelokurovetaL ( 2006 ) distances (the dot- 
ted line in Figure [4| we recover the solid line seen in 
Fig ure [i] 



Newberg et al. (2007) distances are derived using 



photometric parallaxes of main-sequence F-type stars 
with a distance calibration from Sgr BHBs. As above, 
systematic changes in metallicity along the leading arm 
may induce systematic divergence in derived distance. 
Nonetheless, as can be seen from Figure |4] the assump- 
tion of low metallicity produces a degree of concor- 
dance between the two previous determinations and 
the present study. 

3.1 Discussion of substructures 

In addition to the Sgr leading arm, four spatially dis- 
tinct regions of over-density are apparent in Figure [4] 
We identify these as the Virgo Over-density (VOD) 
that has been known from previous studies discussed 
above and new candidate stellar over-densities I term 
the 160° and 180° features, and the Virgo Equatorial 
Stream (VES). 

The first issue we must consider is if these fea- 
tures could result from a background population of 
red clump stars in our sub-giant colour selection win- 
dow, due, for instance, to the Sgr stream. In order to 
produce the Virgo Equatorial Stream, 160° and 180° 
over-densities the Sgr stream red clump stars would 
have to lie at '-^ 110, 65 and ^ 75kpc respectively. 
As evident from numerous studies of Sgr (see above 
discussion, and shown here in Figure |4| the leading 
arm is a factor of two closer at these positions and 
so can not account for the observed over- density via 
its population of red clump stars. In the case of the 
VOD, the distance to its associated RRL over-density 
is in concordance the distances derived here. 

Similarly, to bring about the observed over-densities 
be means of lower main-sequence stars falling in the 
sub-giant colour selection window would require dis- 
tances of approximately 7, 4, and 5kpc for the Virgo 
Equatorial Stream, 160° and 180° over-densities. There 
are no indications from the distribution of RRLs at 
these distances that such over-densities do exist ( |Miceli| 
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[2008|. 

Secondly, the features do appear to be artefacts of 
sampling to a coarse spatial grid. They are not ex- 
punged by a different choice of spatial pixel centroids. 
We have experimented with sub-pixel spatial shifts of 
the grid and recover the over-densities at similar sig- 
nificance. 

The region examined in the present study has been 
the focus of a number of previous studies targeting halo 
substructure. In the following section we discuss our 
results in the context of these previous studies where 
possible. To define the extent of each feature we imple- 
ment a kernel discriminant analysis. The technique is 
a non-parametric smoothing of the spatial distribution 
by an optimally selected multivariate kernel. The re- 
sult is a probability density function for the data that 
retains maximum information content afforded by the 
data (i.e. not over/under smoothed). Such probability 
density function were constructed (using the ks pack- 
age (Duong 2007) in R) for all four features. Contours 
that contain 75% and 50% of feature data (drawn from 
a broad rectangular volume containing the feature) are 
shown projected in two dimensions in Figure [5] The 
corresponding areas of the features are given in Table 

m 

We draw the reader's attention to the fact that 
these are sectional slices through the data cube and so 
as such only show those features visible at a specific 
distance. For this reason the Sgr leading arm does 
not appear as a contiguous element but as a 6-9 de- 
gree wide section in th ese panels. The Monoce r os rin g 
dNewberg et al.||2002| ) (or warp, "Conn et al.l ((20081) 
falls outside the field under consideration here at lower 
galactic latitudes off the right of Figure [4] 

In Figure [6] we show the Hess diagrams for each 
sub- structure population. These diagrams are derived 
by taking the on-source population CMD density and 
subtracting an adjacent off-source CMD density, ap- 
propriately scaled to match the number of stars in 
the on-source field. Here, g and r magnitudes have 
been dereddened according to the reddening maps of 
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Schlegel et al. ( 1998 ) together with a colour-specific 
reddening index ofE(^-r) = 1M2E{B-V) and Ag/E{B- 
= 3.793. The Hess diagrams reveal low density main- 
sequence turnoff features for each sub-structure. For 
example. Figure [g] (top left) shows the go ^ 19 turnoff 
for the VOD. In this example, the more distant Sgr 
leading arm material can be seen in the sub-giants at 
go - 21. 



3.1.1 The Virgo over-density 

The VOD is seen to be a diffuse structure centred at 
RA=194.6°, Dec=+1.8° covering 144 square degrees 
and distances from ^ 12-18 kpc. With a peak signifi- 
cance at a distance of 16 kpc, the projected dimensions 
of the VO D in RA is 5 kpc and 1.7 kpc in De c. We note 
that .Keller et al.| ( [2009t ; [Prior etal] ( |2009t ; [Newberg 
et al.| ( |2007| ) (discussed iDelow) have shown that the 
VOD extends outside the bounds of the SDSS survey 
to the south-west and hence the spatial extent derived 
here is a lower limit. Further more, the significance 
measured along the line of sight for the VOD is also 



Figure 6: Hess diagrams of the population of the 
sub-structures discussed in the present study. We 
have subtracted from the on-source field a popula- 
tion taken from an adjacent off- source field, appro- 
priately scaled to match the total number of stars 
in the on-source field. Lighter greyscale indicates 
population excess. The arrows indicate the ap- 
proximate magnitude of the main-sequence turn- 
off in each case. 



underestimated due to the depth of the structure. This 
-\/)inderestimation occurs because our method for deter- 
mining significance requires the subtraction of a low 
order polynomial fit to the general background. Our 
simulations show that this has the effect of diminish- 
ing the significance of, but not erasing, features that 
exhibit a line of sight depth of p > 0.4 mag. 

The VOD was first detecte d as an over-densi ty in 
RR Lyrae variables (RRLs) by Vivas et al. (2001 ) and 
independently detected as an over-dens ity of F-type 
main-sequence stars by Newberg et al. (2002). Both 
studies describe the centre of the VOD at RA=190°, 
Dec=0° and a heliocentric distance of 18 kpc. From a 
2.3°-w ide band centred on Dec=-1.18°, [Vivas Zinn] 
([2006t find a peak density for the VOD at RA=186° 
and a dist ance of 17kpc. In the study of RRL can- 
didates by Ivezic et al. (2005) a distance of 18 kpc is 
derived for the VOD with a center at RA^ 190° in a 
2.5°-wide band centred on Dec=0°. 

[Duffau et al.[ ([2006[) presents radial v elocities for 
a sample of [Vivas et al. (2001) RRLs and Sirko et aT 
( ,2004| BHBs. They find six of the nine stars in the 
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densest region of the VOD (centred at RA=186° and 
Dec=-1°) share a common radial velocity of +100 kms~^ 
with a velocity dispersion less than measurement un- 
certainties. This moving group was termed by |Duf-| 
|fau et ah] the Virgo Stellar Stream to differentiate it 
from the general region of over-density. Recently, |Vi-| 
vas et al. (2008') finds members of the VOD that share 



We propose that the northern portion of the stellar 
density seen the Juric et al. result (their figure 37) 
is due to the Sgr leading arm rather than the VOD, 
and that this leads to Juric et al. to greatly overesti- 
mate the spatial e xtent of the VOD. In addition, the 
Juric et al.l ( 2008 ) study utilizes photometric parallax 



this common radial velocity (i.e. the VSS) that extend 
as close as 12 kpc. By quantifying the excess in the lu- 
minosity function in VOD fields relative to offset fields. 



Duffau et al. (2006) reports that the VOD covers an 
area in excess of 106 square degrees. 

It was demonstrated by|Newberg et al.|(|2007|) that 
the VOD extends beyond the southern edge of the 
SDSS survey. Newberg et al. shows significant over- 
density that appears a continuation of the VOD in 
one of the th ree SEGUE 'outrigger' scans. |Newberg| 



to determine the distance to the target main-sequence 
stars. As discussed previously, the inherently large 
uncertainties associated with this technique lead to an 
oi;erestimation of the depth of the V OD by [Juric et al.| 
In the case of the Prior et al. ( 2009[ ) study, the spatial 
sampling proved insufficient to resolve the VOD region 
into t he series of distinct components seen by |Keller[ 
et al.|(|2009). 



et al. (2007) examines SDSS spectroscopy in a 1.5° 



The VOD possesses an apparent line-of-sight incli- 
nation on the sky such that the closest region at 12 
kpc is to the eastern extremity (RA'^ 184°) and the 
furthest region at 18 kpc is at the westerly extremity 



radius field centred on the peak of this over-density 
at RA=191° and Dec=-7.8°. They find a marginally 
significant peak in the radial velocity distribution in 
common with the VSS. This suggests that the VSS 
extends across this area. 



(RA= 2 10°). This is in line with the findings of [Keller 
et al.[ (|2009[) who 



Juric et al. (2008) presents a study of the den- 



sity excess of main-sequence F-type stars towards the 
VOD. They present evidence that the VOD covers 
^1000 square degrees centred on RA=192° and Dec=+2° 
and has extensive line-of-sig ht depth (6-20 kpc). The 
study of Prior et al.[ (|2009[) implem ents a technique 
similar to that of Duffau et al. ( 2006| that utilizes the 
excess in the observed luminosity function. In this 
case, the exce ss is relative to Besancon galaxy model 
( [Robin et al.[|2003| ) predictions. This study estimates 
an areal extent of the VOD of 760 square degrees. 



Keller et al. ( 2009 ) matches the technique of Prior 
et al.[ ([2009[)~oveFr X 1° fields from the SDS S DR6 



^ _2009[ ) who, from distances to the RR Lyrae 

members, finds the westerly portion to be at a distance 
of 20 kpc from a centra l region (at RA=194°) of 17 kpc. 
Since An et al. ( [2009 ) finds no significant metallicity 
gradient across the VOD region, this is likely due to a 
true spatial gradient. 

3.1.2 180° over-density 

The 180° over-density (centred at RA=179.8°, Dec=+2.0°) 
is of unresolved depth (< 2 kpc) and covers 45 square 
degrees of sky. At a distances of 19 kpc, the projected 
radius for this feature is 1.5 kpc. The 180 ° over-density 
was first seen in the Keller et al. (2009) study (their 
Feature C) but without significant RRL population a 
distance was not able to be ascribed for the feature. 



Adelman-McCarthy et eT^OOS) and SEKBO ([Keller[ ^-^-^ 1^0° over-density 



m 

^t al. 2008) datasets. The SEKBO data enables an ex- 
tension to the south-west that reveals that the region 
of the VOD encompasses three regions of over-density. 



Keller et al. (2009} Feature A is a linear feature ap- 



proximately 3 wide and 14.5° long centred at (198°, 
-10°) with a distance of 20 kpc (hence linear size of 
1x5 kpc). Feature A is not seen here as it does not 
extend into the SDSS sample. Feature B at (192°, -2°) 
is 3° wide and 10° long, with a heliocentric distance 
of 17 kpc (linear size of 1x3 kpc). It is coincident 
with, and hence identified as, the VOD. A third over- 
density. Feature C, is identified with the 180° feature 
of the present study (discussed below). 

The extent of the VOD seen in the present study. 



and that of Keller et al. 



than that of [Juric et al. 



2009), is s ubstantially smalle r 
2008[) o r [Prior et al.[ ([2009[). 



The study of [Juric et al.[ ( [2008[ ) derives the spatial 
extent of the VOD from the density of stars within 
0.2< g - r <0.3 and 20< r <21, so chosen to se- 
lect main- sequence stars a t a dis tance appropriate for 
the VOD. jNewberg et al. (2007) applies identical se- 
lection criteria, but for bracketing distances (r magni- 
tudes one magnitude brighter and fainter; see [Newberg[ 



et al.| (2007) figures 5 and 6). The results of Newberg 
et al. show that the Sgr stream remains a significant 
contributor to the stellar density over 20 < r < 21. 



The 160° over-density is centred at RA=160.9°, Dec=13.3° 
and has an extent of 36 square degrees. It is located at 
a heliocentic distance of 17 kpc, at which distance the 
approximately circular feature presents a radius of 1.3 
kpc. The 160° over-density is seen projected in front 
of the Sgr leading arm. In Figure[6](top right) this Sgr 
material is apparent as the sub-giant population seen 
at go - 20.2. 

3.1.4 The Virgo Equatorial Stream 

The Virgo Equatorial Stream (VES) is centred at (RA=203.3° 
Dec=+2.7°, r=28 kpc) and covers at least 162 square 
degrees (since it is not clear that it is contained within 
the southern boundary of the SDSS field). It possess a 
small line-of-sight depth of <4 kpc. The corresponding 
projected dimensions of the VES at a mean distance 
of 28 kpc are 8.9x4.4 kpc. 

The VES is situated slightly westward of the VOD 
but lOkpc more distant (A^ = 1.6 magnitudes). Could 
the VES be a metal-rich sub-population of the VOD? 
The above derived distances assume an old (12Gyr), 
metal-poor population ([Fe/H]=-1.7). This is anal- 
ogous to stellar population seen in present-day faint 
dwarf spheroidals in the local neighbourhood. If the 
observed structures formed from a much larger (M > 
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10 Mq) gas-rich system such as Sgr we might expect a 
broad metalhcity distribution function for the debris. 
This is seen in the case of the Andromeda Giant Stel- 
lar Stream believed to ha ve been recently d erived from 
a M - lO^M© satellite ( |lbata et al]|2007 |. However, 
if we assume [Fe/H]=-1.7 for the VOD, the implied 
VES metallicity would be greatly super-solar. In ad- 
dition, the work of An et al . (2009) does not support 
the presence of a substantial spatial metallicity varia- 
tion in this region. Similarly, it is not possible that the 
VES could result from the lower MS of the VOD since 
the VES is only 1.6 mag. more distant that the VOD 
and not the ~ 3 mag. expected for the lower MS. For 
these reasons we conclude that the VOD and VES are 
separate spatial structures. 



4 Luminosities of sub-structure 

The number of sub-giants associated with each of the 
four features can be equated to an integrated absolute 
magnitude for the underlying stellar populations. For 
each pixel we form the sum of the number of stars in 
excess of the local background within our sub-giant 
box over the magnitude range of the feature. An old, 
metal-poor isochron e (as before, 12 Gyr, [Fe/H]=-1.7 
Girardi et al. ( |2004 )) is then populated according to a 
Salpeter initial mass function (with masses> O.O8M0) 
until the observed number of sub-giants is reproduced 
( Keller et al.||2001 ). We then form the integrated ab- 
solute magnitude of this population. Table [l] reports 
the resulting absolute magnitudes for the features. 

The uncertainties in these values encompass the 
stochastic variation in multiple random realisations of 
the model population. There is also some variation in 
these luminosities depending on the assumed metallic- 
ity of the population. If we were to assume a more 
metal-rich population, say [Fe/H]=-0.4 such as seen 
in the core of Sgr, the derived luminosities would be 
0.8 magnitudes fainter. Both the VOD and the Virgo 
Equatorial Stream extend outside the SDSS survey 
area. Consequently the absolute magnitudes we derive 
for these features represent lower bounds. In addition, 
in the case of the VOD we expect that the extent of 
the substructure along the line of sight will result in a 
further underestimation of its luminosity as discussed 
in Sect ion [3. 1.1 1 

For the VO D we derive an absolute magnitude of 
Mg = -8.9. l Juric et aL| (2008) calculates an ab- 
solute magnitude of Mr = —8.0 mag. for the VOD 
from the sum of the excess stellar population between 
18<r<21.5, a mean distance of 10 kpc, and an ex- 
tent of 1000 square degrees. This e quates to a surface 
brightness of 32.5 mag. arcsec"^ 



late an absolute magnitude of My =-11. 9 mag. (area 
of 760 sq uare de grees and mean distance of 19kpc). 
The Prior et al. (2009) study utilizes the luminosity 



As 



Juric et al. (2008) 



state, this is a lower limit to the brightness of the sub 
structure, since it does not include bright giants or the 
population below the magnitude limit (r=21.5). Fur- 
thermore, as we have discussed above, we claim that 



the Juric et al. ( 2008 ) study overestimates the spatial 
extent of the VOD. 

By evaluating the excess in observed luminosity 
functions to the VOD comp ared to the p redictions of 
the Besancon Galaxy model. Prior et al.| ([2009^) calcu- 



function excess from 16.25 <V< 19.75. 

The present study differs in the determination of 
integrated luminosity from these previous studies since 
it does not sum over a magnitude range of main-sequence 
over-density. Rather, the number of excess subgiants 
measured here is the consequence of stellar evolution 
and an IMF thus the luminosities derived here are total 
luminosities for the over-densities in question. 

5 Origins of sub-structures: Sgr 
debris? 

The tenuous nature of the stellar populations of the 
VOD, 180°, 160° and VES sub-structures suggests that 
they have originated from the tidal disruption of one or 
mo re prog enit or galaxies. In the studies of |Johnston et| 
al. ( 2008 ) and Font et al. ( 2008 ) a comparison is made 



between simulated observable properties of a stellar 
halo and the parent galaxy's accretion history. The 
studies conclude the high surface brightness features 
{/iiv < 30 mag arcsec"^) in the outer halo (R>10 kpc) 
usually originate from the accretion of single satellites. 
Furthermore, bright features are almost exclusively de- 
rived from relatively high mass prog enitors (lO'^-lO^ 
M0). Since stellar evolution proceeds longer, and pro- 
cessed material is retained, in more massive galaxies, 
the material accreted is expected to be of moderate 
to high metallicity. In other words, the stellar popu- 
lations of substructures are not expected to resemble 
old and metal-poor dwarf spheroidal satellites nor the 
underlying component of the stellar halo. 

[Johnston et al.| (|2008 ) shows that these expecta- 
tions are broadly consistent with our understanding of 
the outer halo of the Milky Way and M31. In the case 
of the Milky Way we find one surviving satellite of high 
surface brightness (Sgr) and it is composed of mate- 
rial significantly more metal-rich than the background 
halo population. Observations of M31 (see Ibata e t al.| 
(2007) and references therein) reveal a halo dominated 
for R> 20 kpc by the Giant Stellar Stream. The Gi- 
ant Stellar Stream is comprised of relatively metal-rich 
debris and is thought to have been deposited during 
the accretion of a high mass (~ 10^ M0) progenitor 
on an approximately radial orbit. The average surface 
brightness of the Gian t Stellar St re am i s i^v = 30 ±0.5 
mag arcsec"^ ( jlbata et al. 2001 ). |lbata e t al. (2007| 
note six other substructures, some of which possess a 
metallicity distribution function similar to the Giant 
Stellar Stream and are potentially related to the same 
accretion event. Morphologically these substructures 
appear as streams, arcs, shells and irregular structures 
superimposed on a metal-poor halo. 

Given the apparent dominance of Sgr debris in 
outer halo substructure, and the finding from cosmo- 
logical simulations that the probability of additional 
high mass accretion events is low for the Milky Way, 
can the observed substructure of the present study be 
explained by Sgr debris? 

In Figure |7| we compare the best-fitting prolate 
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Table 1: The spatial and luminosity details of the sub-structures described in the present study. In the case 
of the VOD and Virgo Equatorial Stream the sky coverage (and hence the derived absolute luminosity) 
is a lower limit as the features are not enclosed in the SDSS area. 



Object 
Name 


RA 

Center 
(deg.) 


Dec 
Center 
(deg.) 


Sky 
Coverage 
(sq. deg.) 


Surface 
Brightness 
(mag./arcsec^) 


Absolute 
Magnitude 
Mg (mag.) 


Mean 
Distance 
(kpc) 


Virgo 
Over-density 


194.6 


1.8 


144 


30.5 




16 


160° 


160.9 


13.3 


36 


32.2 


-5.6l^oi 


17 


180° 
Virgo 
Equatorial Stream 


179.8 
203.3 


2.0 
2.7 


45 
162 


33.2 
31.7 


_r 1+0.4 
^•-■--0.6 

o O+0.5 


19 
28 



model (£=1.25) from the N-body simulations of Law 
et al.| ( |2005| ) to the over-densities of the present study. 
As we have seen in Section [2.11 the choice of metallic- 
ity does significantly modify the derived distance to the 
material and here we assume a single global metallicity 
of [Fe/H]=-1.7 that is significantly less than the mean 



metallicity of the leading arm material ( jChou et al 



2007[ Fe/H=-l.l). It is also the case that there is a 
strong metallicity gradient along th e leading arm ma - 
terial from the core of [Fe/H]=-0.4 ( |Chou et al.|2007D . 
These caveats acknowledged, some qualitative conclu- 
sions can be made between the distribution of model 
debris particles and the observed sub-giant excess in 
the Sgr stream. 

A clear discrepancy is seen in the model's inabil- 
ity to match the Sgr leading arm: namely that for 
RA<170° the Sgr leading arm material is at a greater 
distance at given RA than the model predicts. The 
stellar over-density may be brought into line with the 
model distribution by assuming a much higher metal- 
licity for the debris ([Fe/H]> —0.4). More problem- 
atically, the metallicity distribution function for the 
debris would have to be increasingly metal-rich as we 
proceed away from Sgr, a conclusion discrep ant with 
existing measurements by Chou et al. (2007|. 



To date no one N-body simulation has been capa- 
ble of satisfact orily reproducin g all of the extant obser- 
vational data (Law et al. 2009|. The angular preces- 
sion of the leading arm indicates a spherical or slightly 



oblate halo ( [Fellhauer et al. 2006| |Martmez-Delgado| 
et al. 2007) whereas, the distance to, and radial ve- 
locity of, the leading arm material indicates a prolate 
halo ( Helmi|2004 ). These conclusions have been drawn 
from simulations that are ax isymmetric. The recent 
triaxial models of Law et al. ([2009) are able to repli- 
cate the data, however as pointed out by |Law et al.| 
([2009 ) these models are somewhat unsatisfactory since 
they are in what is expecte d to be an unstable c onfigu- 
ration. The observ ations of [Prior et al^ ( 2009b| ) mirror 
the conclusions of Law et al.| ( |20 05')~that the current 
models of Sgr disruption are inadequate as they do 
not consider the orbital evolution of Sgr. Furthermore, 
as demonstrated by many simulations (see for exam- 



ple, Bekki ( |2008[ )) it may be necessary to consider the 
evolution of both the luminous and dark matter com- 
ponents of the Sgr dwarf as such treatment modifies 
orbital evolution significantly. 

Clearly additional, and more complex, modelling is 
vital to clarify the Sgr debris streams. Given the large 
discrepancy between the current models and our obser- 
vations it is premature to make any definitive remarks 
about the association of this study's sub-structure to 
the Sgr debris stream. Currently, it is possible that the 
Virgo Equatorial Stream could be associated with the 
younger trailing arm material (see Figure |7| and that 
the VOD, 180° and 164° sub-structures may represent 
over-densities along the older trailing arm. However, it 
must be acknowledged that such possible associations 
may not be upheld by radial velocity determinations 
and revised modelling. 

In addition to critically needed modelling, further 
sky coverage in the southern hemisphere is required to 
trace the spatial extent of these substructures. This 
will soon be provide d by the SkyMap per telescope's 
Southern Sky Survey ( [Keller et al.|2007t . Radial veloc- 
ities and metallicities for samples in the sub-structures 
revealed by this study are needed. They will enable us 
to discern if the features can be accommodated by re- 
vised Sgr models or if they represent the remnants of 
distinct merger events. 



6 Summary 

An examination of the excess of sub-giant stars over a 
subset of the SDSS DR6 reveals a series of extended 
halo over-densities in the survey heliocentric distance 
range of 10 - 40 kpc. The most prominent feature is 
that of the Sagittarius dwarf galaxy (Sgr) leading tidal 
debris arm. We derive a spatial configuration for the 
Sgr leading arm material that is in agreement with 
previous studies. The Virgo Over-density is seen to 
extend from 12-18 kpc, with a mean distance of 16 
kpc, and to cover an area of 144 square degrees. This 
is considerably smaller extent than found in previous 
studies. We contend this is due to our resolution of 
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Figure 7: We compare the Sgr model of Law et al. 
( 2QQ5[ prolate model that best fits the leading arm 
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Figure 3: The excess of sub-giant stars seen in our survey volume. In this figure the excess data cube is 
projected along the declination axis. It shows the prominent leading arm of Sgr tidal debris that enters 
the survey volume at RA=120° and a heliocentric distance of 20 kpc and exits at RA~ 190° at 40 kpc, 
the distance limit of our survey. The four inset panels show the colour-magnitude diagrams of 3° x 3° 
fields that cont ain prominent globul ar clusters. The isochrones shown for M3 and NGC 5466, 5024 and 
5053 are due to Girardi et al. (|2004l) (age=12Gyrs; [Fe/H]=-1.7) placed at the distance modulus derived 
from the peak of sub-giant significance from our technique discussed in Section 2.1 The isochrone for 



Pal3 is shown at a distance modulus of four times that measured by our technique to account for the red 
clump population (see Section 2.1 for details). 
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Figure 4: As in Figure [3l the excess of sub-giant stars seen in our survey volume projected along the 
declination axis. Spatial pixels associated with globular clusters are ignor ed in creating this fig ure (see 
text for details). Overlaid is the locus of the Sgr leading ar m as described by Newberg et al. (2007 dashed 
line with points). The dotted line shows the locus due to 'Bel okurov et al.l (2006) and the solid line is 
the same result assuming a [Fe/H]=-l.l as found by ^Chou et al. ( 2007[ se~text for details). The Virgo 
over-density is seen over RA=184°-210° and distances of 12-18 kpc. Three new candidate over-densities 
are seen, the Virgo Equatorial Stream and the 160° & 180° features. 
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Figure 5: Top: Dashed lines show the distance at which sectional slices through the significance data 
cube are made to reveal the spatial orientation and morphology of the identified substructures. Bottom: 
RA, Dec slices through the excess data cube at the mean distances of each sub-structure. The distance 
range averaged in each panel is given in the top right. The contours around each sub- structure shows 
the extent of the structure. The inner contour contai ns 5 0% of the sub-giant over-density, and the outer 
contour contains 75% of the over-density (see Section 3.1 for details). Slices through the Sgr leading arm 
can be seen in the projected sky significance maps for the VOD, the more distant 160° feature and the 
Virgo Equatorial Stream. 



